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======================================================================

*** Introduction

Energy is a universally conserved quantity normally associated with a
system.

Energy can come in many different forms (kinetic, potential, thermal,
mass, etc), and can be converted from one form to another.

Energy, however, cannot be created or destroyed (hence it is
conserved).

What is energy?

Believe it or not, I've never found a fully satisfactory answer
to that question.  We can define SPECIFIC FORMS of energy very easily
(e.g., K = (1/2)mv^2), but a general definition that encompasses all
forms discovered so far, and yet to be discovered, eludes us.

Energy as an ability to do work?  OK.  But what is work?  A transfer
of energy.  Circular definition.


In this section, we will adopt a somewhat simplistic view that energy
can be associated with individual objects, and in the process we will
treat objects as if they were particles.

In the past we could get away with this, since many of the fundamental
laws, such as Newton's laws and conservation of momentum, can be
applied to systems in the same manner as particles (e.g., F_net = ma
for individual particles; F_ext = Ma_cm for systems).

As we will find, energy is more subtle, and we can run into trouble
when we don't take the system nature of energy into account.

======================================================================

*** Kinetic energy

Recall K = (1/2)mv^2.
* K is a scalar.
* K=0 for objects at rest.
* K>0 for moving objects.
* K increases with increasing |v|,
* K does not depend on direction of v.

SI unit:
  1 J = 1 kg m^2/s^2 = 1 Nm
Other units:
  1 erg = 1 g cm^2/s^2 = 10^-7 J  [cgs unit]
  1 ft lb = 1.356 J  [british unit]
  1 eV = 1.602 x 10^-19 J  [atomic unit]

Note that these units can be used for ALL forms of energy,
not just kinetic energy.

In fact, these two are specifically used in thermodynamics:
  1 cal = 4.187 J
  1 Btu = 1055 J

======================================================================

*** Work-KE theorem [examples: vertical free-fall w/o and w/ air resistance]

Kinetic energy of an object can be changed by exerting a force on
it.  That force is said to do work.  Work done by a force is defined
in such a way that the work-kinetic energy is valid:
  dK = W_net = sum_k W_k
W_k represents F_k's contribution to changes in K.

Example: vertical free-fall without air resistance

[Draw up/down motion of ball thrown and caught from same height.  Label
points 1 (bottom), 2 (top), and 3 (bottom again)]

On the way up 1->2:
  K_1 = max
  K is decreasing from 1->2
  W_g,1->2 < 0
  K_2 = 0

On the way down 2->3:
  K_2 = 0
  K is increasing from 2->3
  W_g,2->3 > 0
  K_3 = max (same as K_1)

Total trip 1->3:
  K_1 = max
  K decreases 1->2 and then increases 2->3
  W_g,1->3 = W_g,1->2 + W_g,2->3 = 0
  K_3 = K_1 -- no overall change

In each case,
  K_f - K_i = W_net = W_g  (only force acting)

Note that work is defined over a time interval, and is additive.

Example: vertical free-fall WITH air resistance

[Same diagram, with pts 1, 2, and 3]

On the way up 1->2:
  K_1 = max
  K is decreasing from 1->2 -- more so than without air resistance
  W_g,1->2 < 0  [calculated from h same way as before -- air resistance non-withstanding]
  W_air,1->2 < 0 as well  [frictional ==> "slowing"]
  W_net,1->2 = W_g,1->2 + W_air,1->2 < 0
  K_2 = 0

On the way down 2->3:
  K_2 = 0
  K is increasing from 2->3 -- less so than without air resistance
  W_g,2->3 > 0
  W_air,2->3 < 0  [frictional -- still "slowing"]
  W_net,2->3 = W_g,2->3 + W_air,2->3 > 0  (gravity wins)
  K_3 = max-ish (< K_1)

Total trip, 1->3:
  K_1 = max
  K decreases 1->2 and then increases 2->3
  W_g,1->3 = W_g,1->2 + W_g,2->3 = 0  (still cancels)
  W_air,1->3 = W_air,1->2 + W_air,2->3 < 0
  K_3 < K_1

Note that W_g is calculated from dh in the same manner, regardless
of other forces that may be acting.  This implies that W_g,1->2 and
W_g,2->3 cancel as they did without air resistance.

What changes as a result of air resistance, is the work done by
air resistance itself -- negative.  The end result is K_3 < K_1,
with
  K_3 - K_1 = W_air,1->3 < 0

======================================================================

*** Work done by constant force, 1-D

Recall: W_F = F dx   (1-D, constant force)
  [show page with derivation; generalize later]


Example: Constant 5N force acting on 5 kg block starting from rest.
What is K_f after a displacement of 1.5m?
[draw diagram]

Solution 1: kinematic
  a = F/m = 1.0 m/s^2 (constant)
  dx = +1.5 m
  v_i = 0
  v_f = sqrt(v_i^2 + 2 a dx)
      = sqrt(0 + 2(1.0)(1.5))
      = sqrt(3) m/s = 1.73 m/s
  K_f = (1/2)mv_f^2 = (1/2)(5.0)(3) = 7.5J

Solution 2: W-KE theorem
  W_F = F dx = (+5.0N)(+1.5N) = 7.5J   [note: F,dx same dir ==> W>0]
  K_f = K_i + W_F = 0 + 7.5J = 7.5J

Note: K = (1/2) mv^2 ==> v = sqrt(2K/m)
  v_f = sqrt(2(7.5)/5) = sqrt(3) = 1.73 m/s

Does K_f depend on m?  Does v_f depend on m?


Example: Constant 5N force acting on 5 kg block again, but not
starting from rest.  What is K_f after displacement of +1.5m?

Ki = 1.0J with vi>0 [diagram]
  W_F = 7.5J (same calculation)
  Kf = Ki + W = 1.0J + 7.5J = 8.5J

Ki = 1.0J with vi<0 [diagram]

Is Kf = 6.5J?  That would suggest that Ki = -1.0J (- for direction).
But K has no direction.  K is a scalar.  And it is never negative.
  Kf = Ki + W = 1.0J + 7.5J = 8.5J

Explain by following actual path of object as 5N force continues to act.

Relate to free fall, with initial velocity zero, down, and up.


Example: Tug-of-war / multiple forces
  start from rest
  F1 = +15N  (constant)
  F2 = -10N  (constant)
  no friction
  dx = +5.0m  (after establishing "who wins")
[diagram]

W1 = F1 dx = +75J   (K_f = +75J?)
W2 = F2 dx = -50J
Wnet = W1 + W2 = +25J
Kf = Ki + Wnet = 0J + 25J = 25J

F1 transfers 75J of energy to block
F2 transfers 50J of energy FROM block
Net transfer is 25J to block

Where does energy come from / go?
Who REALLY won the tug-of-war?

Add friction:
  fk = -2N
==> W_fk = fk dx = -10J   (friction converts 10J to thermal energy)

Wnet = (+75J) + (-50J) + (-10J) = +15J

Block still speeds up, but less so than without friction.

======================================================================

*** Work done by constant force, 3-D
**** Discussion

[diagram with F (constant), dr (from i to f -- draw curved path)
and angle phi between F and dr]

W_F = ?

Draw F_par, F_perp and mention that one of these components
contributes to W_F -- the other does not.

Relate to a vs v:
[diagram with a in direction of F; v in direction of dr]

Draw a_par, a_perp and relate to changes in magnitude and direction
of v.

K = (1/2)mv^2 changes when |v| changes -- direction of v does not
matter.

==> F_par contributes to a_par, and thus to changes in K

W_F = F_par |dr| = |F| |dr| cos(phi) = F dot dr


W_F < 0 if phi > 90deg
W_F > 0 if phi < 90deg
W_F = 0 if phi = 90deg; or if dr=0 (or F=0)
[demo the dr=0 part]

======================================================================

**** Angled force on sliding block example

Block on horizontal ground (no friction) with 120N force, 30.0 deg
below horizontal, initially at rest.
  K_f after 4.50m displacement
[diagram]

W_F = F dot dr = (120N) (4.50m) cos(30 deg) = 468J
Note: F_par = (120N) cos(30 deg) = 104N contributes to a

Are we there yet?  Not quite...

W_N = N dot dr = N dr cos(90deg) = 0   [note: N is not zero, W_N is zero]
W_g = mg dot dr = mg dr cos(90deg) = 0

W_net = 468J + 0 + 0 = 468J
K_f = K_i + W_net = 0 + 468J = 468J

What if there is friction?

W_net = 468J + 0 + 0 + (-...) < 468J
K_f < 468J

Are W_N and W_g always zero?  Let's find out...

======================================================================

**** Work done by constraint forces [normal/tension]

W_N = 0 in previous example.  Is that always true?

Horizontal floor [diagram]:
  W_N = N dr cos(90 deg) = 0

Non-curved incline (angle phi) [diagram]:
  W_N = N dr cos(90 deg) = 0    (make sure you understand 90 deg)

Curved incline [diagram]:

Naive attempt (draw full displacement)
  W_N = N dr cos(phi) = ???
Does not work since N is not constant

Do it right; subdivide the interval and consider one small increment:
  dW_N = N dot dr = 0  (right angle)
  W_N = sum dW_N = 0

This can be done formally using integrals (later).


What about tension?

[pendulum diagram]

Subdivide the interval and consider one small increment:
  dW_T = T dot dr = 0  (right angle)
  W_T = sum dW_T = 0


Bead on a curved wire...

Same thing...


Is W_N or W_T ever NOT zero?

(1) Surface moves (N) / pivot point move (T)
N or T vs dr is no longer necessarily at a right angle

(2) Compression (N) or stretch (T) motion cannot be neglected
(collisions / rope suddenly goes taut)

======================================================================

**** Work done by gravity

[diagram: block(i) with mg
  --> 3 block(f)'s, one angled up, one angled down,
      and one horizontal -- W_g < 0, > 0, = 0 indicated]
Relate work sign to angle between dr and mg

W_g = mg dot dr
    = -mg y-hat dot (dx x-hat + dh y-hat)
    = -mg dh
Note that horizontal part of displacement is not relevant.

======================================================================

**** Example: free-fall and frictionless slide

Example: m=5.0kg, dh=-20m, vi=0 (free-fall) [diagram], vf=?

Kinematics:
  vf^2 = vi^2 + 2(-g)dy
  ==> vf = sqrt(-2gdh) = 20 m/s (actually -20m/s)

W-KE theorem:
  W_g = -mg dh = +980 J
  K_f = K_i + W_g = 980 J
  v_f = sqrt(2 K_f/m) = 20 m/s

Note: W-KE does not provide direction, but must be down


Example: m=5.0kg, dh=-20m, vi=0 (frictionless slide) [diagram], vf=?

W-KE theorem:
  W_g = -mg dh = 980 J, as before
  W_N = 0  (argued earlier)
  K_f = K_i + W_net = 980 J
  v_f = 20 m/s (to the right this time)

Note: same final speed as free-fall

What if there is friction?  Then K_f = 980 J - |W_fr|

Kinematics:
  a = F_net/m = (N + mg)/m = ??
  N needs to be calculated -- it is not constant
  requires knowledge of incline shape (parametric)
  seems to require a differential equation
  quite difficult

If the question were dt = ? (instead of v_f=?), quite difficult

If the slide involved friction (mu_k given), problem would be
VERY difficult:
  W_N is still zero, but |N| needed to calculate |f_k| (not constant)
  requires knowledge of incline shape

======================================================================

**** Work done by lift force

Lifting an object up against gravity requires positive work.
Consider the following example:

[diagram: 3.0 kg block lifted vertically 15m by tension
 include force diagram]
  W_T = ?

Naive attempt:
  W_T = + T dh = (mg)dh = 441 J

Correct answer, but faulty reasoning.  T is not constant and
cannot be definitively evaluated.  It might be reasonable
to expect that T = mg "on average", but what does that mean?

Correct attempt; backdoor solution using W-KE theorem:
  W_g = -mg dh = -441 J  (true since mg is constant -- never mind T)
  dK = K_f - K_i = 0-0 = 0  (assuming start/end at rest)
  dK = W_net = W_g + W_T
  0 = -441 J + (W_T)
  W_T = +441 J

When work done can be calculated for all forces except one, and
motion is known (dK can be calculated), W-KE theorem can be used
to calculate the unknown work done.


This trick is often done with air resistance and/or friction on a
curved slide.

Previous example:
  m=5.0kg, dh=-20m, vi=0, vf=18m/s (instead of 20m), W_fk = ?

  K_f = (1/2)mv_f^2 = 810 J  (instead of 980 J)
  dK = W_g + W_N + W_fk
  810 J = 980 J + 0 + W_fk
  W_fk = -170 J
Who cares if fk is not constant.

======================================================================

**** Work done by (constant) friction

1-D: [diagram shows block with displacement dr and friction fk]
  W_fk = fk dot dr = fk d cos(180 deg) = -fk d  (d = |dr|)

What about curved path?
[diagram showing top view curved path]

Naive attempt: (indicate overall dr and fk at start)
  W_fk = fk dot dr = ??
Does not work: fk (vector) is not constant, even if |fk| is.

Correct attempt:

Subdivide the time interval.  At each step,
  dW = fk dot dr = fk |dr| cos(180 deg) = -fk |dr|

Now sum over all steps.  Since |fk| is constant (even though
direction changes), we can factor it out and get
  W = -fk sum(|dr|) = -fk d
where d = DISTANCE along path (not |dr|!!!)

If |fk| is not constant, integrals must be used.

======================================================================

**** Example [slide down frictionless incline, then stop by friction]

Diagram setup:
  Block (5 kg) starts at top at rest (20m curved incline)
  curved incline is frictionless
  at bottom, long straight patch with friction (muk = 0.40)

How far does block slide along bottom before coming to rest?

Slide down the incline part (already solved):
  vi = 0 (Ki = 0)
  W_g = -mg dh = +980 J
  W_N = 0
  Kf = Ki + W_net = 980 J = (1/2) m vf^2
==> vf = sqrt(2Kf/m) = sqrt(392) m/s =approx= 20m/s

Slide along the bottom / Newton's 2nd law solution:

[force diagram: N, mg, fk]
  N = mg
  fk = ma = muk N = muk mg
==> |a| = muk g = (0.40)(9.8 m/s^2) = 3.92 m/s^2
  vf^2 = vi^2 + 2a dx
  dx = [vf^2 - vi^2]/2a
     = [0^2 - sqrt(392)^2]/2(3.92)
     = 50 m  (exact)

Slide along the bottom / W-KE:
  fk = muk mg = 19.6 N  (still need N2L for |N|)
  W_fk = -fk d
  W_g = 0  (horizontal)
  W_N = 0  (isn't it always?)
  dK = W_fk + W_g + W_N
  -980 J = W_fr + 0 + 0
  -980 J = -(19.6 N) d
  d = 50 m

Note, algebraically, that d = h/muk.

======================================================================

*** Work done by variable force
**** General work formula (derivation)

[diagram: show block following i->f curved path -- show different
forces at different places]

How do we calculate W_F if F is not constant?

Subdivide the path into small pieces, and calculate work done
for each small segment -- then add results.

If segments are small enough, F can be treated as constant over the
small segment.

For one segment:
  dW = F dot dr   [indicate in diagram]

Total work is obtained by integration:
  W = int_path F dot dr

Note that if the force is constant, then
  W = int_path F dot dr
    = F dot int_path dr
    = F dot Dr

If the path is along a coordinate axis (ex: x-axis, from x_i->x_f),
then dr = dx x-hat and
  W = int_path F dot dr
    = int(xi->xf) F dot dx x-hat
    = int(xi->xf) F_x dx
This can be evaluated if F_x is known as a function of x.

======================================================================

**** Parametric paths

In general, a path can be defined parametrically: r(alpha)
  [diagram: show path with pts labeled by alpha and draw r(alpha)]

The parameter is a coordinate which determines location along the
path.  Parameterization is arbitrary.

Consider a small segment described by the parameter interval
[alpha,alpha+dalpha]:
  [diagram: path, r(alpha), r(alpha+dalpha), dr]

In this case,
  dr = r(alpha+dalpha) - r(alpha)
     = [r(alpha+dalpha) - r(alpha)]/dalpha  dalpha
     = [dr/dalpha] dalpha

This allows us to convert the path integral into an integral
over alpha:
  W = int_path F dot dr
    = int(alpha_i->alpha_f) F dot (dr/dalpha) dalpha

This integral is easily done if F is known as a function of alpha,
which will surely be the case if F is known as a function of r:
  F = F(r(alpha))

Back to x-axis:
  r(x) = x x-hat  [x_i,x_f]  -->  dr/dx = x-hat
==> dr = (dr/dx) dx = (x-hat) dx = dx x-hat


If r(t) is known, path can be parameterized by alpha=t, and thus
  W = int(ti->tf) F dot (dr/dt) dt
    = int(ti->tf) (F dot v) dt
This would allow us to calculate work as long as F is known as a
function of time, which would be the case if it is known as a function
of r, v, and/or time.

Time parameterization may not be convenient if motion is not known
ahead of time.

======================================================================

**** Work done by spring (1-D)

[picture of block attached to spring / wall at x=xi (indicate Fs)
 immediately below that, picture at x=xf (indicate larger Fs)]

Work done by the spring as the block is moved from xi->xf?

Hooke's Law: F(x) = -kx x-hat

W = int(xi->xf) F_x(x) dx
  = int(xi->xf) (-kx) dx
  = -(1/2)kx^2 |_xi->xf
  = -(1/2)k[xf^2 - xi^2]

Note that this result is correct, even if x<0 for some or all of
the path.

What if xf < xi [draw situation]?

Then dr = |dx| (-x-hat), but dx<0, so dr = dx x-hat (without sign)
and we get the same result.

Note that the work done by a force that ONLY depends on position will
be opposite if path is reversed.  This is true in general for any path
(3-D).

======================================================================

**** Example: Work done along a straight line [coordinate-directed / diagonal]

]draw diagrams throughout]

Consider:
  F(x,y) = A(x+y) x-hat + Bxy^2 y-hat

(0,0)->(c,0):

Use x as parameter:
  r(x) = x x-hat + 0 y-hat = (x,0)
  dr = (dr/dx) dx = dx x-hat

  F(r(x)) = F(x,0) = Ax x-hat + 0 y-hat

  W = int(0->c) F(x,0) dot dx x-hat
    = int(0->c) F_x(x,0) dx
    = int(0->c) Ax dx
    = (1/2)Ax^2 |_(0->c)
    = (1/2)Ac^2

(c,0)->(c,d):

Use y as parameter:
  r(y) = c x-hat + y y-hat = (c,y)
  dr = (dr/dy) dy = dy y-hat

  F(r(y)) = F(c,y) = A(c+y) x-hat + Bcy^2 y-hat

  W = int(0->d) F(c,y) dot dy y-hat
    = int(0->d) F_y(c,y) dy
    = int(0->d) Bcy^2 dy
    = (1/3)Bcy^3 |_0->d
    = (1/3)Bcd^3

(0,0)->(c,0)->(c,d):
  W = (1/2)Ac^2 + (1/3)Bcd^3


(0,0)->(0,d):

Use y as parameter:
  r(y) = 0 x-hat + y y-hat = (0,y)
  dr = dy y-hat

  F(r(y)) = F(0,y) = Ay x-hat + 0 y-hat
  W = int(0->d) F(0,y) dot dy y-hat
    = int(0->d) F_y(0,y) dy
    = int(0->d) 0 dy
    = 0    [note F is perpendicular to dr at all points]

(0,d)->(c,d):

Use x as parameter:
  r(x) = x x-hat + d y-hat = (x,d)
  dr = dx x-hat

  F(r(x)) = F(x,d) = A(x+d) x-hat + Bxd^2 y-hat

  W = int(0->c) F(x,d) dot dx x-hat
    = int(0->c) F_x(x,d) dx
    = int(0->c) A(x+d) dx
    = A((1/2)x^2 + d*x) |_(0->c)
    = A(c^2/2 + cd)

(0,0)->(0,d)->(c,d):
  W = 0 + A(c^2/2 + cd)
Note that this is different from (0,0)->(c,0)->(c,d)


(0,0)->(c,d) direct line:

Use alpha parameter:
  r(alpha) = c*alpha x-hat + d*alpha y-hat   alpha:0->1
  dr = (dr/dalpha) dalpha = (c x-hat + d y-hat) dalpha
Note the direction dr is along the line joining (0,0 and (c,d).

  F(r(alpha)) = F(c*alpha,d*alpha)
              = A(c+d)alpha x-hat + Bcd^2 alpha^3 y-hat

  W = int(0->1) F(r(alpha)) dot dr
    = int(0->1) [A(c+d)alpha x-hat + Bcd^2 alpha^3 y-hat] dot [(c x-hat + d y-hat) dalpha]
    = int(0->1) [A(c+d)alpha c  +  Bcd^2 alpha^3 d] dalpha
    = Ac(c+d) (1/2)alpha^2  +  Bcd^3 (1/4)alpha^4  |_0->1
    = (1/2)Ac(c+d) + (1/4)Bcd^3

======================================================================

**** Example: Work done along a circular path

Consider again:
  F(x,y) = A(x+y) x-hat + Bxy^2 y-hat

(+R,0)->(-R,0) along semi-circle (Quadrants I and II):

Use phi as parameter:
  r(phi) = R(cos(phi) x-hat + sin(phi) y-hat)   phi:0->pi
  dr = (dr/dphi) dphi
     = R(-sin(phi) x-hat + cos(phi) y-hat) dphi  [= R phi-hat dphi]

  F(r(phi)) = F(Rcos(phi),Rsin(phi))
            = AR(cos(phi)+sin(phi)) x-hat
                + BR^3 cos(phi) sin(phi)^2 y-hat

  W = int_(0->pi) F dot dr
    = int_(0->pi) [AR(cos(phi)+sin(phi)) x-hat + BR^3 cos(phi) sin(phi)^2 y-hat]
                     dot R(-sin(phi) x-hat + cos(phi) y-hat) dphi
    = int_(0->pi) [-AR^2(cos(phi)+sin(phi))sin(phi)  +
                    BR^4 cos(phi)^2 sin(phi)^2] dphi
    = -AR^2[0 + pi/2] + BR^4 [pi/8]
    = -(pi/2)AR^2 + (pi/8)BR^4

Note:
  int(0->pi) sin(phi)cos(phi) dphi = 0
  int(0->pi) sin(phi)^2 dphi = pi/2
  int(0->pi) sin(phi)^2 cos(phi)^2 dphi = pi/8


Contrast with straight line (+R,0)->(-R,0):

  W = int(R->-R) F_x(x,0) dx   [note: dx<0]
    = int(R->-R) A(x+0) dx
    = (1/2)x^2 |_R->-R
    = 0

======================================================================

**** Work done by a central force (F(r) r-hat) (radial/tangential/general path)

A central force can be expressed as a function of position in
the following form:
  F(r) = f(r) r-hat
  ^^^^   ^^^^
 vector scalar

where f(r) is the (signed) radial component of F, and can only
depend on the radial coordinate.  A central force obeys spherical
symmetry around the origin (force center).

The inverse square forces associated with universal gravitation and
Coulomb force between charges are special cases of this.


Work done by a central force along a radial path [diagram: ri->rf]:

Parameterize the path with the radial coordinate:
  r(r) = r r-hat  r:ri->rf
  dr = (dr/dr) dr = dr r-hat

  W = int(ri->rf) [f(r) r-hat] dot [dr r-hat]
    = int(ri->rf) f(r) dr

For an inverse square force (f(r)=K/r^2), the result is
  W = K(-1/r) |_ri->rf = -K(1/rf - 1/ri)


Work done by a central force along a circular path [diagram: phi_i->phi_f]:

Parameterize the path with the angular coordinate:
  r(phi) = R (cos(phi) x-hat + sin(phi) y-hat)
  dr = dr/dphi dphi = R (-sin(phi) x-hat + cos(phi) y-hat) dphi
     = R phi-hat dphi

  W = int(phi_i->phi_f) [f(R) r-hat] dot [R phi-hat dphi]
    = int(phi_i->phi_f) 0 dphi    [F and dr are at right angles]
    = 0

This result can be generalized for any motion contained entirely
within a sphere of radius R around the origin.
  [try to draw this]
At all points along the path, F and dr will be at right angles, and
W=0.


Work done along a general path:   [draw this in 2-D / polar]
  r(alpha) = r [sin(theta)[cos(phi) x-hat + sin(phi) y-hat] + cos(theta) z-hat]
           = r r-hat    [r, theta, phi are all functions of alpha]
  dr = (dr/dalpha) dalpha
     = [r' r-hat + r theta' theta-hat + r theta phi' phi-hat] dalpha
where r'=dr/dalpha, theta'=dtheta/dalpha, and phi'=dphi/dalpha

This is often written
  dr = dr r-hat + (r dtheta) theta-hat + (r theta dphi) phi-hat
where
  dr = r' dalpha; dtheta = theta' dalpha; dphi = phi' dalpha

  W = int(alpha_i->alpha_f) F dot dr
    = int(alpha_i->alpha_f) [f(r) r-hat] dot [r' r-hat + (...)theta-hat + (...)phi-hat]
    = int(alpha_i->alpha_f) f(r) r' dalpha  [tangential parts = 0]
    = int(ri->rf) f(r) dr    [substitution: alpha->r]

======================================================================

**** Revisiting work done by constraint and friction forces

We already addressed these cases informally.  Now let's consider the
formal derivations involving integrals.  In each case, it proves to be
most convenient to use time parameterization.

[draw diagram]

For an object moving along a surface, N and v are at right angles.
It follows that
  W_N = int_path N dot dr
      = int(ti->tf) N dot v dt    [time parameterization]
      = int(ti->tf) 0 dt
      = 0

Note: this is NOT true if the surface moves, or if compression motion
cannot be neglected -- in this case N and v might not be at right
angles.

The same thing happens with tension (with same caveats) [draw diagram].

Kinetic friction fk along a 2-D path (assume constant |fk|):
  W_fk = int_path fk dot dr
       = int(ti->tf) fk dot v dt
       = int(ti->tf) -|fk| |v| dt    [fk and v are at 180 deg angle]
       = -|fk| int(ti->tf) |v|dt
       = -|fk|(total distance)

======================================================================

**** Proof of Work-KE theorem

Forget we ever heard of W-KE theorem, and define the work done by
a force by the integral:
  W_F = int_path F dot dr
      = int(ti->tf) F(t) dot v dt   [time parameterization]

Can we prove the W-KE theorem?  [Of course, otherwise why are we here?]

  K = (1/2)mv^2 = (1/2) m (v dot v)
  dK/dt = d[(1/2)m v dot v]/dt
        = (1/2)m 2v dot dv/dt
        = ma dot v
        = F_net dot v
        = sum_k F_k dot v

Then
  DK = K_f - K_i
     = int(ti->tf) (dK/dt)dt
     = int(ti->tf) (sum_k F_k dot v) dt
     = sum_k int(ti->tf) F_k dot v dt
     = sum_k W_k
     = W_net

QED

======================================================================

*** Power

Power is defined as the rate of work done.

Average power:
  P_avg = (work done) / (total time)

If the force is constant over the time interval in question, then
  P_avg = W/Dt = (F dot Dr)/Dt = F dot v_avg

Instantaneous power:
  P = lim(dt->0) P_avg = F dot v
This is true regardless of whether F is constant or not.

Note that
  dK/dt = sum_k F_k dot v = sum_k P_k = P_net
This is the work-KE theorem, expressed in terms of rates, rather
than overall changes during a finite time interval.

Note that a given force will deliver more power if the object upon
which it acts is moving faster (along the direction of the force).
This makes sense, since work is being done faster as more distance
is covered.

SI units:
  1 W = 1 J/s

Other units:
  1 hp = 746 W   [roughly the power output of one horse]
  1 ft lb/s = 1.356 W

Any energy transfer rate can be expressed in terms of power.  For
example, rate of heating an object can be expressed in J/s or cal/s
or Btu/s, etc.

A 60W lightbulb uses 60J of electric energy per second.

======================================================================

*** Potential energy
**** Gravity and spring demonstration

[demo 220a]
Lift object against gravity -- must do positive work -- and yet
no change in KE.  Where did lift energy go?

[demo 220b]
Hold one end of spring while stretching spring on other end --
must do positive work -- yet no change in KE.  Where did energy go?

A: Energy supplied by my force is converted into potential energy.

When I allow a return to the original configuration, KE comes back.
[complete 220a and 220b videos]

In both cases, amount of KE that comes back is equal to the work
done initially.

Energy is not lost -- it is converted into potential energy.

Potential energy always comes back when objects are returned to their
original positions.

======================================================================

**** Conservative forces

Examine gravity case in detail:

[diagram: A below B, object moved from A to B; W_g,A->B < 0]

As I lift object against gravity, gravity also acts and does
negative work -- explains apparent energy disappearance.

When object returned from B->A [indicate in diagram], gravity does
positive work: W_g,B->A, thus restoring the KE that is apparently
lost during the lift phase.

What makes this work is the fact that
  W_g,B->A = -(W_g,A->B)

This is true regardless of
(1) points A and B
(2) path followed
(3) time at which it happens

A conservative force is a force for which
  W_B->A = -(W_A->B)
for all points A and B, paths, and timing.

Equivalent definition of conservative force:
  W_A->A = 0  for any closed path   [prove equivalence]

Gravity and spring forces are conservative.

[repeat demos 220a and 220b]

======================================================================

**** Potential energy

For conservative forces, PE defined as a function of position, as follows:
  dU = -W
  U_B - U_A = - W_A->B
  W < 0 ==> dU > 0  (U increases as force does negative work)
  W > 0 ==> dU < 0  (U decreases as force does positive work)

Example: gravity

[diagram: A->B, m=2.0kg, dh=1.0m]

  W_g,A->B = -mgdh = -19.6J
  dU = U_g,B - U_g,A = -W_g,A->B = +19.6J
  U_g,A = 0  (reference point)
  U_g,B = 19.6J

  W_g,B->A = -mgdh = +19.6J = -W_g,A->B
  dU = U_g,A - U_g,B = -19.6J  [consistent]

Note importance of W_B->A = -(W_A->B) condition.  Without it,
U cannot be defined consistently.

Analogy:
  energy = money
  gravity = bank
  W = transaction
  KE = cash
  PE = bank account balance

What if U_g,A=5J?  U_g,A=-3J?  Reference values do not matter -- as
long as differences are correct.

Note that if U_A = 0, then U_B for any point B will equal W_B->A.

Potential energy = future work
  (potential for future work when object is returned to its U=0
   reference position)

======================================================================

**** Non-conservative forces

Are all forces conservative?

Consider friction.
[demo 220c]

[diagram: object moved from A->B; W_fr,A->B < 0, W_fr,B->A < 0 again!]

I do positive work to overcome friction A->B.  Where does energy go?

Does it come back?

[complete demo 220c]

No.  More positive work done -- energy disappears again.

Analogy continued:
  friction = thief/mugger

There is no "frictional potential energy".

Energy is converted into a different form -- mostly thermal energy.

======================================================================

**** Gravitational potential energy

To calculate potential energy from force:
(1) Pick two arbitrary points (one point can be fixed)
(2) Calculate work done by force moving from one point to the other.
    Choice of path should not matter (conservative force).
(3) Calculate DU = -W.
(4) Choose reference point and place one of the two points at that point.
(5) Solve for the potential energy at the other (arbitrary) point.


Consider terrestrial gravity (F_g = mg):

[diagram: m travels from h_i to h_f along some path]

  W_g = -mg Dh = -mg(hf - hi)
  DU_g = -W_g = mg(hf - hi)

Choose reference height: U_g=0 when h=0 (simplest choice)

  U_g(h) - U_g(0) = mg(h - 0)
  U_g(h) = mgh

Notes:
(1) U_g=0 when h=0 [based on choice of reference height]
(2) U_g increases as h increases.
    higher height ==> more potential to do positive work
(3) mg itself points from higher potential energy to lower potential energy.
(4) U_g can go negative, if h<0.  This represents a "negative account
    balance" (a "loan") -- object needs to be lifted to reach U_g=0J.

Example revisited:
  m = 2.0 kg
  g = 9.8 m/s^2
  hi = 0.0 m  (ref)
  hf = 1.0 m

  U_gi = mg hi = 0 J
  U_gf = mg hf = 19.6 J
  DU = +19.6J = -W_g
  etc...

======================================================================

**** Elastic (spring) potential energy

[diagram: block connected to spring connected to wall --
 block is moved from 0 to x (Fs points left)]

We will choose reference point (U_s=0) to be at equilibrium (x=0).
Choose one point to be 0, other point is arbitrary (x):

  W_0->x = int(xi->xf F_x(x) dx   [path along x-axis]
         = -(1/2)k(x^2 - 0)       [calculation done earlier -- show page]
         = -(1/2)kx^2
  U_s(x) - U_s(0) = -W_0->x = +(1/2)kx^2
  U_s(x) = (1/2)kx^2

Notes:
(1) This formula is valid for x <= 0 (compression) as well as x >= 0
(2) Potential energy is minimum at equilibrium -- spring does negative
work whenever it is stretched or compressed.
(3) U_s >= 0 always -- consequence of choice of reference point
(4) F_s = -kx x-hat points from higher potential energy to lower
potential energy, as it always does.


Example:
  Fs = 4.0N when x=10cm ==> k=40N/m
  Us = (1/2)kx^2 = 0.20 J

When stretching spring from 0->10cm, applied force must do +0.20 J
of work to offset spring force (doing -0.20 J of work).

When x=20cm,
  Us = (1/2)kx^2 = 0.80 J   (x *2 ==> U_s *4)

Note that 0.60 J of work is needed to stretch spring by another
10cm -- force got stronger.

======================================================================

**** Potential energy for inverse square forces

Central force:
  F(r) = f(r) r-hat
  W_(ri->rf) = int(ri->rf) f(r) dr   [see "Work done by a central force"]

Note: This is independent of path -- central forces are conservative.

Inverse square force:
  f(r) = K/r^2
  W_(ri->rf) = -K(1/rf - 1/ri)

Where to place reference point (U=0)?  At r=0?  Doesn't work.

Set reference point at r=infinity (1/r=0):
  U(r) - U(inf) = -W_(inf->r)
                = +K(1/r - 1/inf)
  U(r) = K/r

Notes:
(1) When K>0 (repulsion), U>0 and decreases as r->infinity.
(2) When K<0 (attraction), U<0 and increases (becomes less negative)
    as r->infinity
(3) In either case, force points from higher PE -> lower PE

Gravity (extra-terrestrial):
  F_g = -G m1 m2/r^2 r-hat  [Newton's universal law of gravitation]
  ==> U_g(r) = -G m1 m2 / r  (negative and increases as r increases)

Electric:
  F_e = k q1 q2/r^2 r-hat  [Coulomb's law]
  ==> U_e(r) = k q1 q2 / r  [sign depends on charges]

======================================================================

**** Status of constraint (normal/tension) and exerted (push/pull) forces

We already know that gravity (both terrestrial and universal) and
elastic (spring) forces are conservative, and that frictional forces
are not conservative.

What about normal and tension forces?  This would include push/pull
forces.  Are these forces conservative?

Given the limitations of object view and our current discussion of
conservative vs non-conservative forces, there is no easy answer to
this question.

Issues:

1. Normal and tension forces often do zero work.  In such cases, it
almost doesn't matter whether they are regarded as conservative or
not.

2. Normal and tension forces are often modelled as compression/stretch
forces, which are "spring-like".  In this context, it makes sense to
regard them as conservative.

3. Sometimes, however, there is a non-conservative nature to such
compression/stretch forces.  This can occur in extreme circumstances,
such as (inelastic) collisions.  Same thing can happen with (regular)
springs if they are over-stretched or over-compressed.

4. When there is a sentient being (or otherwise complex system) connected
with this force (e.g., push/pull), the force may be selectively applied
or withdrawn.

Case in point: focus on my lift force in demo 220a.

5. This last point actually demonstrates one of the limitations of our
current "object-oriented" view of energy, as the same issue can befall
gravity if the earth is moved (consider demo 220a again, but imagine
that the earth is moved away before the eraser returns to its original
position).

When we consider a force on an object to be conservative or not
conservative, we are assuming that the "environment" that exerts the
force is not changing.  When we consider system view, we will address
this.


Considering the issues above, the best course of action (at least
for now) is to assume that normal/tension forces are not conservative,
and to not assign a potential energy associated with those forces.

======================================================================

*** Conservation of mechanical energy
**** Vertical free-fall discussion [w/o and w/ air resistance]

Consider vertical free-fall.  [diagram]

On the way up: K decreases as U_g increases
On the way down: K increases as U_g decreases

U_g is defined in such a way that changes in K due to gravity are
offset by changes in U_g.

It therefore follows that, in the absence of air resistance,
  E_mech = K + U_g = (1/2)mv^2 + mgh   [definition of mechanical energy]
does not change.

Mechanical energy, defined above, is conserved during free-fall
processes.

In this case, K is converted to U_g as the object rises; U_g is
converted back to K as the object falls; but the total does not
change.

What if we take air resistance into account?

In that case, K will decrease on the way up, partially due to gravity
(accounted for by increases in U_g) AND partially due to air
resistance (not accounted for in U_g at all).  It follows that E_mech
will overall decrease by the work done by air resistance.

On the way down, K increases, but not by as much as U_g decreases.
Again, the difference is accounted for by W_nc.

It turns out that mechanical energy
  E_mech = K + U
is conserved in any situation where the only forces acting are
conservative.

If non-conservative forces act, then mechanical energy is not
conserved, but changes according to
  DE_mech = W_nc

The changes in mechanical energy represent conversion of mechanical
energy to/from other forms of energy.

Frictional forces tend to convert mechanical energy into (mostly)
thermal energy.

======================================================================

**** Derivation of dE_mech = W_nc

Consider object acted on by one force that is conservative.

[diagram: object moving from i->f -- one force acting]

  DK = W_net = W_F   [W-KE theorem]
  DU = -W_F          [definition of U]
==>
  DE_mech = DK + DU = (W_F) + (-W_F) = 0
==>
  E_mech,i = E_mech,f
  K_i + U_i = K_f + U_f

Mechanical energy is conserved.


Now suppose object is acted on by several forces, all conservative.

  DK = W_net = sum_k W_k   [W-KE theorem]
  W_k = -DU_k   [definition of U_k -- one term for each force]
  U_total = sum_k U_k   [total PE]
  DU_total = sum_k DU_k = sum_k (-W_k) = -W_net
  E_mech = K + U_total = K + sum_k U_k   [def of mechanical energy]
  DE_mech = DK + DU_total
          = W_net - W_net
          = 0

Mechanical energy is conserved in this case as well.


Now suppose object is acted on by several forces, some conservative
and some non-conservative.

  DK = W_net = W_cons + W_nc = -DU_total + W_nc
  DE_mech = DK + DU_total = W_nc

Mechanical may not be conserved, but changes by W_nc.


What if non-conservative forces act, but none of them do any work?

Then W_nc = 0 and mechanical energy will still be conserved.

======================================================================

**** Example: Box launched from cliff

Setup: [diagram]
  10 kg box
  launched from top of cliff at undetermined angle
  height of cliff = 40 m

  U_top = mgh = 3920 J  (ref height at bottom)
  U_bot = 0 J
  dU = -3920 J
  W_g = +3920 J  (confirmed by W_g = -mg dh)

Suppose v_i = 20 m/s
  K_top = 2000 J
  E_mech,top = 5920 J  (same value throughout motion; assume no air resistance)
  E_top = E_bot ==> K_bot = 5920 J  (confirmed with dK = W_net)

Suppose W_air = -400 J (air resistance)
  dE_mech = W_nc = -400 J
  E_mech,bot = 5520 J
  K_bot = 5520 J  (confirmed with dK = W_net)

======================================================================

**** Example: Pendulum swing

Pendulum swing:  [diagram]
  L = 4.0 m
  m = 2.0 kg
  v_bot = 8.0 m/s
  ignore air resistance

[demo 230]

Is E_mech conserved?   Yes (W_nc = 0).

Q: at phi = 60 deg, what is v?

* K_i + U_i = K_f + U_f
* K based on speed
* U based on position
* i: position/speed known; both K_i and U_i can be calculated
* f: position known: U_f can be calculated
* solve for K_f ==> v_f can be calculated

  K_i = 64 J
  U_i = 0 J  (ref height)

[Derive h = L(1-cos(phi)]

  h_f = 2.0 m  (at phi=60 deg)
  U_f = 39.2 J
  K_f = 24.8 J
  v_f = 5.0 m/s

Same speed every time phi=60 deg is reached on either side

General consequence of E_mech cons: same speed every time same
position reached.


Q: phi=? when v=0

  K_i = 64 J; U_i = 0 J  (no change)
  K_f = 0 J
  U_f = 64 J
  h_f = 3.27 m
  phi_f = 79 deg

Won't reach higher at any time.

Note: U <= E_mech always, since K >= 0.


What if air resistance is taken into account?

E_mech slowly decreases over time:
* speed at 60 deg is < 5.0 m/s (if it gets there) and decreases each time.
* phi_max < 79 deg and decreases each time

General feature of oscillations.

Note the irony: definition of E_mech includes U, which involves
conservative forces, but conservative forces have no effect on
E_mech (precisely *because* changes in K due to such forces are
offset by changes in U).  However, non-conservative forces can
affect E_mech.

======================================================================

**** Example: Spring gun (horizontal)

Spring gun launch (horizontal):  [diagram]
  m = 12 g
  x_i = -3.2 cm
  k = 750 N/m  (7.5 N/cm)
  v_f = ? (after launch)

  K_i + U_s,i = K_f + U_s,f

  K_i = 0 J
  U_s,i = 0.384 J
  U_s,f = 0 J
Note: U_s actually depends on spring configuration -- not mass location
  K_f = 0.384 J
  v_f = 8.0 m/s  (don't shoot this at anyone; even a toy gun can injure)

Effects of friction / air resistance.

Effects of massive spring.  [demo 240]

======================================================================

**** Example: Vertical spring launch

Setups: [side-by-side diagrams]
  initial spring, nothing on it
  m=8.00kg placed on spring -- compresses 10.0 cm
  compress an additional 30.0 cm and release
  how high does it reach before coming to rest again?

[demo 240]

  mg = 78.4 N
  k = 784 N/m
  K_i = U_g,i = 0 J  (ref height)
  U_s,i = 62.7 J  (note: x_i = 40.0 cm)
  K_f = U_s,f = 0 J
  U_g,f = 62.7 J
  h_f = 80.0 cm  (40.0cm above uncompressed height)

Details:
* From -40.0cm to -10.0cm (|F_s| > |F_g|):
    conversion: U_s -> K, U_g
* From -10.0cm to 0.0cm (|F_s| < |F_g|):
    conversion: K, U_s -> U_g
* From 0.0cm to +40.0cm (F_s = 0)
    conversion: K -> U_g

Note: dK = 0 overall, even though K changes during interval

Notice how complex this process is.  But a conservation principle
is still a conservation principle: U_g,f = U_s,i when all is said
and done because mechanical energy is conserved during this process.

Air resistance?  Massive spring?

What if calculated h_f is below the uncompressed height?

======================================================================

**** Example: Spring-mass oscillations

Diagram: horizontal spring-mass oscillations

[see "Spring-mass oscillations" in forces.txt]

Mass released from rest from x=A.
Solution:
  x(t) = A cos(omega*t); omega = sqrt(k/m) = angular frequency
  v(t) = dx/dt = -A*omega sin(omega*t)

Is mechanical energy conserved?

U_s = (1/2)kx^2
    = (1/2)k [A cos(omega t)]^2
    = (1/2)kA^2 cos(omega t)^2
K = (1/2)mv^2
  = (1/2)m [-A omega sin(omega t)]^2
  = (1/2)m omega^2 A^2 sin(omega t)^2

omega = sqrt(k/m) --> m omega^2 = k

E_mech = K + U_s
       = (1/2) kA^2 [sin(omega t)^2 + cos(omega t)^2]
       = (1/2) kA^2

Time dependence disappears -- E_mech is constant (conserved)

Note that E_mech transfers from potential (at x = +/- A; v=0) to
kinetic (at x=0, v = v_max = +/-A*omega) and back again as
oscillations continue.

Note that conservation of mechanical energy is directly related
to the fact that the emplitude of the oscillation remains
constant.  We will see later that the presence of friction will
cause A to decrease over time [demo 450].


What about vertical oscillations?

[See "Mass hanging on spring (equilibrium)" in forces.txt]

[diagram: mass hanging from vertical spring]

x = x0 + x'
x = position (+ = down) relative to original equilibrium (at 0)
x0 = mg/k = adjusted equilibrium position
x' = position relative to adjusted equilibrium

U_tot = U_s + U_g
      = (1/2)kx^2 + mgh
      = (1/2)kx^2 - mgx  (ref height at original equilibrium)
      = (1/2)kx^2 - k x0 x
      = (1/2)k[x^2 - 2 x0 x]
      = (1/2)k[(x-x0)^2 - x0^2]
      = (1/2)k x'^2 - (1/2)k x0^2

First term = potential energy of spring as if its new equilibrium were the
actual spring equilibrium

Second term = constant adjustment, which has no physical consequences

This is consistent with the idea, introduced earlier, that
spring+gravity is equivalent to spring with an adjusted equilibrium
length.

Second term = -(1/2)k x0^2 = mg[(-1/2)x0]

Readjusting reference height downwards by (1/2)x0, so that
  h = (1/2)x0 - x
will get rid of this second term.

Indeed,
U_tot = (1/2)kx^2 + mgh
      = (1/2)kx^2 + mg[(1/2)x0 - x]
      = (1/2)kx^2 + k x0[(1/2)x0 - x]
      = (1/2)k[x^2 - 2 x0 x + x0^2]
      = (1/2)k(x-x0)^2
      = (1/2)k x'^2

The full solution for these spring-mass oscillations are
  x(t)-x0 = x'(t) = A cos(omega t)
  v(t) = -A omega sin(omega t)

The same calculation reveals that mechanical energy (which includes
spring and gravity PE) is conserved.

Note: The result for vertical spring-mass oscillations, and the
calculation U_s+U_g above, can be used to solve the previous problem
in the event that the mass is not launched from the vertical spring.

======================================================================

**** Escape velocity

***** Maximum height for terrestrial gravity

Object (mass m) launched upwards (initial velocity vi).

Maximum height (h_max) reached?

Ignoring air resistance and using terrestrial gravity:
  K_i + U_gi = K_f + U_gf
  (1/2) m vi^2 + 0 = 0 + mg h_max   [ref height at launch]
  h_max = vi^2 / (2g)
It seems that what goes up must eventually come down.

See previous calculation using free-fall:
  vi = 5.0m/s ==> h_max = 1.28m
Works well in this case.

What about this case?
  vi = 10^4 m/s ==> h_max = 5.1x10^6 m = 5100km

(1) Ignoring air resistance would be unreasonable here, but put that
aside.

(2) R_e = 6.38x10^6 m.  Gravity becomes markedly weaker over this
distance range.

***** Maximum height with universal gravity

Continue to ignore air resistance, but now take variation in g into
account.
  F_g = -GMm/r^2 r-hat
where
  G = 6.673x10^-11 N m^2/kg^2
  M = 5.98x10^24 kg  [earth's mass]
  R = 6.38x10^6 m  [value of r at earth's surface]
==> U_g = -GMm/r  [ref point at infinity]

  g0 = GM/R^2 ==> GM = R^2 g0  [g0 = 9.80m/s^2 = g at earth's surface]
  U_g = -m g0 R^2/r

  K_i + U_gi = K_f + U_gf
  (1/2) m vi^2 - m g0 R^2/R = 0 - m g0 R^2/(R+h_max)
  (1/2) m vi^2 = m g0 R^2 [1/R - 1/(R+h_max)]
               = m g0 R [1 - 1/(1+h_max/R)]
Define,
  gamma = vi^2 / (2 g0 R)
        = [vi^2/(2 g0)] / R
          ^^^^^^^^^^^^^----- value of h_max assuming terrestrial gravity

Then,
  gamma = 1 - 1/(1+h_max/R)
  1 + h_max/R = [1-gamma]^-1
  h_max = R [1/(1-gamma) - 1]

If vi = 10^4 m/s, then
  gamma = vi^2 / (2 g0 R)
        = 0.7997
  h_max = R [1/(1-gamma) - 1]
        = 3.992 R
        = 25500 km   (much higher)

***** Escape velocity

What happens when gamma = 1?  This occurs when
  vi = sqrt(2 g0 R) = 11.2 km/s
At this point h_max -> infinity and the object never returns.

This speed is known as the escape velocity.

Note that
  U_gi = -GMm/R = -m g0 R
  K_i = (1/2) m vi^2
and so
  gamma = -(K_i / U_gi)

If gamma = 1, then U_gi = -K_i and therefore
  E_mech = K_i + U_gi = 0
As the object gets farther away from the earth, U_g will
increase (i.e., become less negative) towards zero, and K->0.

The kinetic energy will not reach zero until r->infinity.

If gamma > 1 (vi > escape velocity), then
  E_mech = K_i + U_gi > 0
and the kinetic energy will approach a finite (positive) value as
r->infinity.  There will be kinetic energy left over when the
object gets very far away from the earth.

By contrast, if gamma < 1 (vi < escape velocity), then
  E_mech = K_i + U_gi < 0
and the kinetic energy will reduce to zero at a finite value of
r (where U_g is still negative).  The object will never reach
infinity, and will instead come back down.

So, what comes up usually comes down, but may not if vi > 11.2km/s


Note: If gamma << 1, then
  1/(1-gamma) =approx= 1 + gamma   [first-order expansion]
and so
  h_max = R[1/(1-gamma) - 1]
        = R[1 + gamma - 1]
        = R gamma
        = vi^2 / (2g0)
as calculated initially.

======================================================================

**** Example: Lift example revisited

[diagram: 3.0 kg block lifted vertically 15m by tension
 include force diagram]
  W_T = ?

Recall original solution:   [show previous page]
  dK = W_g + W_T
  0-0 = -mg dh + W_T
  0 = -441 J + W_T
  W_T = +441 J

Alternative solution:
  dK = 0-0 = 0
  dU_g = mg dH = +441 J
  W_T = W_nc = dE_mech = dK + dU_g = +441 J

Work done by tension force is transferred to block's potential energy

======================================================================

**** Example: Slide down incline and stop revisited

Diagram setup:
  Block (5 kg) starts at top at rest (20m curved incline)
  curved incline is frictionless
  at bottom, long straight patch with friction (muk = 0.40)
How far does block slide along bottom before coming to rest?

Indicate in diagram (calculated from before):
  At top: K = 0, U_g = 980 J  (E_mech = 980 J)
  At bottom: K = 980 J, U_g = 0  (E_mech = 980 J)
  After stop: K = 0, U_g = 0  (E_mech = 0)
  fk = muk mg = 19.6 N

Kinematics solution (bot->stop):
  vi = sqrt(2 K_bot/m) = sqrt(392) m/s
  vf = 0
  a = fk/m = -3.92 m/s^2
  d = (vf^2 - vi^2)/2a = 50 m

W-KE solution (bot->stop):
  dK = W_net = W_g + W_N + W_fk
  -980 J = 0 + 0 + (-fk d)
  d = 980/19.6 = 50 m

E_mech solution (bot->stop):
  dE_mech = W_nc = W_fk
  -980 J = -fk d
  d = 50 m

Alternative solution (top->stop):
  dE_mech = W_nc = W_fk
  -980 J = -fk d  (friction only acts along the bottom)
  d = 50 m

======================================================================

*** More on potential energy
**** Force as the derivative of potential energy (1-D)

Situation:
  Object can move in 1-D (x-axis)
  U(x) is given
  What is F?

Suppose object moves from x -> x+dx (diagram).
  W_F = F dx = -dU = -[U(x+dx) - U(x)]
  F = -[U(x+dx) - U(x)]/dx = -dU/dx

Important consequence:

No matter what v or t is, whenever object is at x, the force
acting on it must be F(x) = -dU/dx, as calculated from the potential
energy function U(x).

Conservative forces must be position-only dependent.

Friction, in all of its forms, is disqualified, since it is invariably
velocity dependent (even sliding friction, whose direction depends on
the direction of v).

Another consequence:

Adding a constant value to U(x) does not change F.  Changing the
reference point / reference value has no physical effect.

Examples:

U_s(x) = (1/2) kx^2
==> F_s(x) = -dU_s/dx = -kx

U_g(y) = mgy   [+y is up]
F_g(y) = -dU_g/dy = -mg

======================================================================

**** Potential energy graphs (1-D)

A lot can be learned about the motion of an object acted on by a
conservative force from a graph of U(x).


[example hilly graph]

F(x) = -dU/dx
(1) steeper graph --> stronger force
(2) force points from high U to low U  [indicate force directions]

Stationary points of U --> F=0 (equilibrium points)

Local maxima of U --> unstable equilibrium
Local minima of U --> stable equilibrium


U_g(h) = mgh -- proportional to h

Strong analogy with motion of object sliding (no friction) along a
hill shaped like the U(x) graph

Not a complete analogy, since real motion along a hilly terrain
involves vertical and horizontal motion, so the distances are
different.


Now draw E_mech as a horizontal line.

  K = E_mech - U
so speed of object at x is related to difference between E_mech line and
U(x).

Allowed region of travel: all points x where U(x) <= E_mech

When U(x) = E_mech, K=0 and object stops.  If dU/dx is not zero there,
it must have approached from below (U<E_mech side), and will turn
around and go back to that side.

[Explain back-and-forth motion in diagram]


If friction is present, E_mech will decrease over time -- restricting
motion.  Eventually, object will be confined to one "potential well",
and oscillate back and forth, slowly decreasing in amplitude until the
object comes to rest "at the bottom" (stable equilibrium).

======================================================================

**** Force as the gradient of potential energy (2,3-D)

(x,y,z) -> (x+dx,y,z)

W_F = F dot dx x-hat = Fx dx
    = -[U(x+dx,y,z) - U(x,y,z)]
Fx = -[U(x+dx,y,z) - U(x,y,z)]/dx = -dU/dx [y,z=constant]  -- partial derivative

(x,y,z) -> (x,y+dy,z)

W_F = F dot dy y-hat = Fy dy
    = -[U(x,y+dy,z) - U(x,y,z)]
Fy = -[U(x,y+dy,z) - U(x,y,z)]/dy = -dU/dy [x,z=constant]

(x,y,z) -> (x,y,z+dz)

Fz = -dU/dz  [x,y=constant]

Gradient:
  grad(U) = (dU/dx) x-hat + (dU/dy) y-hat + (dU/dz) z-hat

--> F = -grad(U)

Examples:
  U_g(x,y,z) = mgy  (+y is up)
  F_g = -[(dU_g/dx) x-hat + (dU_g/dy) y-hat + (dU_g/dz) z-hat]
      = -mg y-hat

Radial potentials:
  U(x,y,z) = U(r)
  r = [x^2 + y^2 + z^2]^(1/2)  [radial coordinate]

  dr/dx = (1/2)[x^2 + y^2 + z^2]^(-1/2) (2x)
          = x/r
          = (r-hat)_x
  Fx = -dU/dx
     = -dU/dr dr/dx
     = -U'(r) (r-hat)_x
Ditto for y,z.

  grad(r) = r-hat
  F = -grad(U) = -(dU/dr) grad(r) = -U'(r) r-hat

U(r) = A/r --> F = A/r^2 r-hat

======================================================================

**** Potential energy contour plots (2,3-D)

Just like 1-D, a lot can be learned about the motion of an object
acted on by a conservative force from a graph of U(x,y) in 2-D, or
U(x,y,z) in 3-D.

An explicit graph of U(x,y) requires a 3-D plot (x, y, and U axes),
and can be difficult to create by hand, although computer graphing
programs exist that can plot it for us:
  [Use grapher program to generate 3-D plot: z = 3/(1 + x*x + y*y)]

An explicit graph of U(x,y,z) requires a 4-D plot, which is especially
difficult to draw, or even imagine.

Fortunately, graphical representations exist for U(x,y) and U(x,y,z)
that are easier to draw and visualize: contour plots.


May be familiar with contour maps that indicate elevation using contour
lines:
  [draw a simple, but non-symmetric "hill" -- 10m, 20m, 30m, 40m, 50m]

* Each line (curve) represents points at the same elevation.

* Steepness can be inferred from closeness of contour lines

Contour plot of U(x,y) (or U(x,y,z)) is similar
  [draw similar example: 10m->10J, etc]

* Each contour represents points with the same PE value.

* U(x,y) contours are curves; U(x,y,z) contours are surfaces.

* grad(U) points "uphill" -- F = -grad(U) points "downhill".

* F (and grad(U)) is always perpendicular to contours --
F dot dr must be zero for dr along the contour, as no work can be
done by F if U remains constant.

* Steepness is related to |F| -- force is larger in regions where
contours are closer together.


Explicit examples:

Terrestrial gravity (contour surfaces are horizontal planes -- equally
spaced).

Inverse-square (contour surfaces are spheres centered at force center
-- not equally spaced).  Discuss escape velocity motion.


Example (pre-drawn) involving peaks and wells.

Motion can be deduced based on PE graph -- imagine a ball rolling
over the terrain.

E_mech defines a limiting value of U -- determines allowed region.

Important difference with 1-D: Edge of allowed region may not be
reached because of the possibility of lateral motion (object can
turn around without stopping).

If friction acts, E_mech slowly decreases, narrowing the allowed
region.  Eventually, object will be confined to one well and will
come to rest at the bottom.

======================================================================

**** Conservative forces revisited
***** Four equivalent conditions for conservative forces

A force is conservative if and only if it is position-only dependent
and satisfies one of the following equivalent conditions:

(1) Work done by force depends only on initial and final point, and
not on the path followed by the object:
  int(ri->rf,path 1) F dot dr = int(ri->rf,path 2) F dot dr  [diagram]

(2) Work done by force around a closed loop is always zero
  closed-int F dot dr = 0   [diagram]

(3) Force can be written as (minus) the gradient of a potential energy
function:
  F = -grad(U)

(4) The curl of F is zero everywhere
  curl(F) = 0

We already know about the first 3 conditions, specifically that (1) and (2)
are equivalent and each imply (3).  The fact that (3) implies (1) (or (2))
can be derived with the help of vector calculus (later).

We will be spending much of our time in these upcoming videos on condition (4).

Condition (4) appears to come from nowhere, but the curl (defined later)
can be directly related to integrals around small closed loops, as we
will show.

It has the advantage that it can be easily checked.

======================================================================

***** Conservative forces are curl-free

Suppose F is a conservative force (satisfies (2)).

Consider a small (ccw) rectangular loop aligned with the x and y axes:
dx x dy centered at (x,y,z) [diagram].

The loop consists of four segments, assumed small enough that we can
approximate the work done along each segment with a single term:

  W_F = F(x,y-dy/2,z) dot [+dx x-hat]      seg 1
        + F(x+dx/2,y,z) dot [+dy y-hat]    seg 2
        + F(x,y+dy/2,z) dot [-dx x-hat]    seg 3
        + F(x-dx/2,y,z) dot [-dy y-hat]    seg 4
      = - dx [Fx(x,y+dy/2,z) - Fx(x,y-dy/2,z)]   segs 1&3
        + dy [Fy(x+dx/2,y,z) - Fy(x-dx/2,y,z)]   segs 2&4
      = dx dy {
         [Fy(x+dx/2,y,z) - Fy(x-dx/2,y,z)]/dx
         - [Fx(x,y+dy/2,z) - Fx(x,y-dy/2,z)]/dy }
      = dx dy [ dFy/dx - dFx/dy ]
      = dx dy [ grad_x F_y - grad_y F_x ]   <-- this looks like a cross product
      = dx dy [grad x F]_z
               ^^^^^^^^----- curl of F, z-component

If F is conservative, then W_F=0, which requires
  (grad x F)_z = 0

If we shift x,y,z -> y,z,x and consider a rectangle aligned
with the y and z axes, we can do a similar calculation and show
that
  W_F = dy dz [ dFz/dy - dFy/dz ] = dy dz [grad x F]_x = 0

Considering a rectangle aligned with the z and x axes, we can
also conclude
  W_F = dz dx [ dFx/dz - dFz/dx ] = dz dx [grad x F]_y = 0

Putting everything together leads to
  curl(F) = grad x F = 0

The fact that curl(F) = 0 is a sufficient condition for a force to be
conservative (i.e., (4) implies (2)) can be derived with the help of
vector calculus (Stoke's Theorem).  We will address this later.

======================================================================

***** 2-D force example revisited

Is condition (4) always satisfied?

For 1-D forces (e.g., aligned along the x-axis, so that Fx depends
only on x and Fy=Fz=0), then yes, since none of the cross derivatives
are non-zero.  1-D position-only dependent forces are always
conservative.

For general (2,3-D) forces, the answer is no.  We used this example
from before (see "Work done along a straight line").
  F(x,y) = A(x+y) x-hat + Bxy^2 y-hat

We got different answers for total work along
  (0,0)->(c,0)->(c,d) and
  (0,0)->(0,d)->(c,d)
This force is non-conservative.

The curl (only z-component is non-zero) is given by
  curl(F)_z = dFy/dx - dFx/dy
            = d(Bxy^2)/dx - d(A(x+y))/dy
            = By^2 - A
It is not zero everywhere.

======================================================================

***** Sufficiency of condition (3)

In this video, we will demonstrate that condition (3) implies
condition (1).

To that end, suppose that
  F = -grad(psi)
where psi(r) is some (scalar) function.

At this point, we have not established that F is conservative, and
psi should not be identified with PE at this time.

Consider a path described parametrically by r(alpha), where alpha
varies from alpha_i to alpha_f.  We wish to show that
  W_F = int F dot dr
depends only on the end points r_i=r(alpha_i) and r_f=r(alpha_f),
and not the path, thereby establishing (1).

(d/dalpha)psi(r(alpha))
  = (d/dalpha) psi(x(alpha),y(alpha),z(alpha))
  = dpsi/dx dx/dalpha + dpsi/dy dy/dalpha + dpsi/dz dz/dalpha
  = (dpsi/dx,dpsi/dy,dpsi/dz) dot (dx/dalpha,dy/dalpha,dz/dalpha)
  = grad(psi) dot dr/dalpha

W_F = int F dot dr
    = int(alpha_i->alpha_f) F dot dr/dalpha dalpha
    = int(alpha_i->alpha_f) [-grad(psi)] dot dr/dalpha dalpha
    = -int(alpha_i->alpha_f) (d/dalpha)psi(r(alpha)) dalpha
    = - psi(r(alpha))|_alpha_i->alpha_f
    = -[psi(r(alpha_f)) - psi(r(alpha_i))]
    = -[psi_f - psi_i]

Indeed, W_F is independent of path, and so the force is conservative.
Furthermore, U(r)=psi(r) is (one possible) potential energy function
associated with this force.

[see "Potential energy for inverse square forces",
"Work done by a central force", and "Force as the gradient...",
for inverse square force example]

======================================================================

***** Sufficiency of condition (4)

In this video, we will demonstrate that condition (4) implies
condition (2), thereby completing the proof that all four conditions
are equivalent.

The key here is to first recognize that curl(F) gives us information
about closed integrals of F around SMALL loops centered at
individual points [see "Conservative forces are curl-free"].

If curl(F)=0 at all points (condition (4)), then all such integrals
will be zero.

To prove (2), we need to extend this result to large loops.

Consider the following two adjacent rectangular loops:
  [diagram: both loops ccw arrows]

If we integrate both loops in the same direction, we find that
  int(big loop) F dot dr = int(loop 1) F dot dr + int(loop 2) F dot dr

If we add up the integrals along each side of each smaller loop, the
integrals along the common edge connecting the two loops cancel.  The
remaining left-over edges add up to the larger loop.

This additive property of closed-loop integrals can be extended to any
number of smaller loops combined to form one larger loop.

Consider an arbitrary (large) loop.  Subdivide into many very small
loops.  [diagram of large loop rectangularized, with ccw arrows]

  int(big) F dot dr = sum_k int(loop k) F dot dr
                    = sum_k 0   (since curl(F)=0)
                    = 0
establishing (2).

This proof can stated more rigorously with the help of Stoke's
Theorem (from vector calculus):
  int(closed-loop) F dot dr = int curl(F) dot da  (da is a vector)
where the RHS is integrated over the area bounded by the closed loop.

curl(F) dot da represents the closed-loop integrals along one of the
small rectangles.

da is made into a vector in order to pick off the correct component
of curl(F) [see "Conservative forces are curl-free"].

The entire RHS will be zero if curl(F)=0 everywhere, thereby
establishing (2).

======================================================================
NOT-YET

*** Advanced situations
**** Time down a frictionless curved incline

parametric path given: r(alpha): 0 -> alpha_f

The trick is to determine |v(alpha)| using W-KE.
  W_g = -mg dh = mg(y(0) - y(alpha))
  K_f = K_i + W_g = (1/2) mv_i^2 + mg(y(0) - y(alpha))
  v_f = sqrt(2K_f/m) = sqrt(v_i^2 + 2g(y(0) - y(alpha)))
Note that m cancels out -- mass does not affect motion

Now consider interval [alpha,alpha+dalpha]:
  |v| = sqrt(v_i^2 + 2g(y(0) - y(alpha)))
  |dr| = |dr/dalpha| dalpha
       = sqrt(x'(alpha)^2 + y'(alpha)^2) dalpha
  dt = |dr|/|v|
     = sqrt(x'(alpha)^2 + y'(alpha)^2) dalpha / sqrt(v_i^2 + 2g(y(0) - y(alpha)))

To get the total time down the incline, we simply integrate the above
expression:
  T = int(0->alpha_f) [sqrt(x'^2 + y'^2) / sqrt(v_i^2 + 2g(y(0)-y))] dalpha
This can be done numerically.

Spreadsheet example...

======================================================================
NOT-YET

**** Curved incline with friction

Preliminaries:
  r(alpha) given.
  r' = dr/dalpha; r'' = d(r')/dalpha
  beta = |r'| = sqrt(r' dot r')
  (beta^2)' = 2 beta beta' = (r' dot r')' = 2 r'' dot r'
  ==> beta' = (r'' dot r')/beta

  tau = unit tangent vector = r'/beta
  kappa = curvature vector
        = tau'/beta
        = [r'/beta]'/beta
        = [r''/beta - r'/beta^2 beta']/beta
        = [r''/beta - r'/beta^3 r'' dot r']/beta
        = [r'' - r'(r'' dot r')/beta^2]/beta^2
        = [r'' - tau(r'' dot tau)]/beta^2

...

======================================================================

*** Limitations of object view

So far we have been studying energy and energy conservation from the
point of view that energy is associated with an object (i.e., "object
view").

It turns out that object view has limitations.


(1) External conservative forces...

When considering a conservative force acting on an object, we
tacitly assumed that the object moves in some static "force field".
--- Whatever exerts that force is omnipresent, and never changes.

E.g., with gravity, the earth remains in place to exert the force.

Suppose we lift an object A->B against gravity, remove the earth,
and then move the object back B->A.

In this case
  W_g,A->B = -mg dh
  W_g,B->A = 0  [earth is gone!!]

No longer appears to be conservative.

The work done on the earth in moving it away needs to be
accounted for.

It also needs to be noted that U_g depends on the position of the
object RELATIVE TO THE EARTH -- not its absolute position in space.

System view will clarify this.


Consider a free-fall situation with an object dropped from rest.  If
we measure K_obj in a downward moving reference frame, K will
initially decrease as object falls.  U_g also decreases (since
distance to earth decreases), and so
  E_mech = K + U_g
appears not to be conserved.

The KE of the earth needs to be accounted for.

System view will clarify this as well.


(2) Treatment of objects as if they were simple particles...

Kinematics:
  r(t) = ___; v = dr/dt; a = dv/dt
  r_cm(t) = ___; v_cm = dr_cm/dt; a_cm = dv_cm/dt   OK

Newton:
  F_net = ma
  F_ext = M a_cm    OK

Momentum:
  p = mv
  P_sys = M v_cm   OK

Energy:
  K = (1/2) mv^2; W = int F dot dr
  K_sys =?= (1/2) M v_cm^2; W_F =?= int F dot dr_cm   NOT TRUE!!!!!

A system cannot be simply reduced to its center-of-mass
motion for calculations of K and W.  Relative motion contributes
to KE and can affect work calculations.

Energy can hide!!!

System and atomic view will clarify this.

Energy associated with small atomic motions is the subject of
Thermodynamics.

Thermodynamics: The price we pay for the right to summarize...

======================================================================
